Depolarizing field and "real" hysteresis loops in nanometer-scale ferroelectric films 
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We give detailed analysis of the effect of depolarizing field in nanometer-size ferroelectric capacitors 
studied by Kim et al. [Phys. Rev. Lett. 95, 237602 (2005)]. We calculate a critical thickness of 
the homogeneous state and its stability with respect to domain formation for strained thin films of 
BaTiOa on SrRuOa/SrTiOs substrate within the Landau theory. While the former (2.5nm) is the 
same as given by ab-initio calculations, the actual critical thickness is set by the domains at 1.6nm. 
There is a large Merz's activation field for polarization relaxation. Remarkably, the results show 
a negative slope of the "actual" hysteresis loops, a hallmark of the domain structures in ideal thin 
films with imperfect screening. 

PACS numbers: 77.80.Dj; 77.80.Fm; 85.50.Gk 



Noh et al. 0, 0, have recently published a series of 
seminal experimental studies of very thin, down to 5nm, 
ferroelectric (FE) BaTiC>3 capacitors. They obtained 
hysteresis loops at frequencies up to 100kHz and stud- 
ied retention of a single domain (SD) state in various 
external fields. They have applied high external field to 
obtain polarization saturated state, perhaps a SD one, 
and observed its relaxation when the external field Eq 
were reduced below a certain value Eq = Eq t . Noh et 
al. have approximately identified this field as a depolar- 
izing field due to incomplete screening by electrodes and 
claimed that it coincides with the one estimated from 
electrostatics. We show below that such an interpreta- 
tion does not apply and present a consistent interpreta- 
tion, which reveals important features of domain struc- 
ture in electroded thin FE films: (i) we find deviations 
from the Merz's empirical formula[^| for thickness and 
exposition time dependence of the activation field for do- 
main wall movement and/or domain nucleation in very 
thin films; (ii) replotting the hysteresis loops as function 
of a field in the ferroelectric reveals a negative susceptibil- 
ity of multidomain films governed by the electrostatics, 
predicted some time ago |5j; (hi) finally, we argue that 
domain structure in the thinnest ferroelectric films with 
electrodes may be different from that in thick films. 

In their study, Kim et al.\j^ followed the idea by Mehta 
et ai.|f| that the incomplete screening of the ferroelectric 
bound charge by electrodes leads to a depolarizing field 
inside the ferroelectric that promotes domain nucleation 
and movement and limits the polarization retention. Kim 
et al. speculated that the depolarizing film can be iden- 
tified with an external field Eq = Eo r that stops the 
polarization relaxation. They claimed that E$ r is very 
close to the depolarizing field Ed- This is incorrect, how- 
ever, since Ed was estimated under finite external field 
Eo from formula by Mehta et al. who did only the case 
of short-circuited capacitor (Eo = 0). If they were right, 
the Merz activation field for the domain motion and/or 
nucleation in their ultrathin samples would have been 
negligible, which is not the case (see below). 



It is easy to find the external field Eq = Eob that com- 
pletely compensates the depolarization field (i.e. corre- 
sponds to zero field in FE). This is the point where the 
field in FE changes sign from negative to positive with 
regards to the polarization. The homogeneous field in 
the monodomain ferroelectric is 

E f = (E -2PX/e e e l)/(l + 2X/e e l). (1) 
= Eq + Ed, 

where Eo = U/l is the external and 

E d = -PLo/eol, L = 2A/e s (2) 

the depolarizing field, Lq the characteristic length scale 
in electrodes, and we have used the fact that in a metal 
the screening length is small, X/e e l <C 1. The relation 
between the polarization and the external field Eo is then, 
using the equation of state Q renormalized by strain : 



A h P + BP 3 + CP 5 + FP 7 = 



Eo 



l + 2X/e e V 



A 



2A 



e {e e l + 2A) 



A 



2A 

tQtel' 



(3) 



(4) 



Here we have taken into account that the first coefficient 
may be renormalized by the additional boundary condi- 
tions (ABC), A = A+(2a + (3)/l .10], but in the present 
case the ABC effect is not important, since the variation 
of the spontaneous polarization on thickness, where the 
(2a + j3)/l term enters, is small, see Fig. la (inset). The 
equation of state for spontaneous polarization P s takes 
the form 



A + BP 2 + CP* + FP° = 0, 



(5) 



which is readily solved analytically. Since we know 
from the data the polarization at zero external field, 
Po = P(E = 0)0, we can indeed easily find P s 
from We see from Eq. JIJ that Ef — when 
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FIG. 1: (a) The external field E ob (where E f = 0) and 
Eor where the relaxation of polarization starts in 5 nm thick 
filmQ. Inset: the spontaneous polarization P s and the ex- 
trapolated Pq 1]. (b) The measured -P(-Bo) and the "actual" 
P(Ef) hysteresis loops. 



Eab = 2AP S / (eo£ e l), the external field for multidomain- 
single domain boundary. 

Finite depolarizing field does not imply that the do- 
mains will be detected after a certain observation time. 
Indeed, — 1450kV/cm in 5 nm sample, which is 
about 60% larger than the relaxation field E$ r identi- 
fied by Kim et al. This means that during the observa- 
tion time t re iax = 10 3 s the domains begin to form only 
when there is a field Ef = — (490 ± 70) kV/cm opposite 
to the polarization (at Eq = E$ r — 910kV/cm for 5nm 
sample, raw data from Kim et a/.|l(). The activation field 
strongly depends on FE film thickness for long applica- 
tion time (~ 10 3 s). Indeed, in the I = 30 nm FE film the 
activation field is practically zero, Fig. la. It is interesting 
to see how the activation field depends on the application 
time. To this end, we have replotted the hysteresis loops 
taken at 2 kHz as a function of a field in the ferroelectric 
E f , P = P{E f ), Figs.lb,2 Note that we used Eq. © to 
calculate Ef, which applies to homogeneous state. One, 
however, can apply it to an MD part of the loops around 
P ~ too, since the domains are narrow and the field 
in the bulk is approximately homogeneous. Even in 5nm 
film the domain width a = 2.2nm (see below), so this 
approximation should be semiquantitative. The use of 
Eq. J2I is justified for finding the activation field for do- 
mains (the ends of a horizontal parts of the loop). We see 
that the activation field of similar magnitude is observed 
at all thicknesses. This field is about the same for 5 nm 
sample as observed for much longer 10 3 s of application 
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FIG. 2: The measured P(E ) and the "actual" P(E f ) hys- 
teresis loops: (a) film thickness I = 9 nm, (b) I =30 nm. 

time. 

Comparing the data by Kim et al with Merz Q, one 
expects that the switching time is close to the application 
time mentioned above. We rewrite the Merz's empirical 
formula in the form 

E a = { a \nr)/l, (6) 

where r is the switching time and a the coefficient. It 
describes a weak dependence of the activation field on the 
application time for I = 5nm as well as strong reduction 
of the field with a sample thickness at a large application 
time rather well. However, weak dependence of E a on 
the film thickness in the case of small application times 
is not captured by ©• 

We find using Eqs. that the electric field 

in short-circuited 5nm sample is Et(Eo = 0) = 
— 1200 kV/cm, exceeding the magnitude of the estimated 
activation field. This means that in a short-circuited 
sample SD state relaxes quicker than in 10 3 s. If the 
value of the activation field is defined by the thickness 
only and not by properties of electrodes or an electrode- 
film interface, one can speculate about the properties of 
electrodes, which can facilitate a smaller field in a short- 
circuited sample and a longer, at least 10 3 s, retention of 
a SD state. We have found that for d = 5nm such an 
electrode should have A/e e < 0.033A. Since in pj this 
value is about 0.1 A, it does not seem impossible to find 
such an electrode. Alternatively, thinner films may show 
longer retention, as Eq. © suggests. 

An unusual specific feature of replotted loops is that 
they all have a negative slope, most pronounced at I = 
5nm. The exact value is affected by error bars, but 
there is a strong indication that it is characteristic of 
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all samples. The negative slope has been predicted some 
time ago for an ideal ferroelectric plate between perfect 
metallic electrodes with a voltage drop across thin dead 
layers 0: it is a hallmark of domain structure governed 
mainly by electrostatics. There are no dead layers in the 
present films Q, but the voltage drop happens across a 
screening layer with the thickness A in an electrode, with 
an identical result and apparently beyond the precision 
of Ref. Q ■ Comparing theory with the data, the expres- 
sion for the dielectric constant, given by Eq. (31) of Ref. 
PJ, can be simplified to 

e/ « -eJ/X, (7) 

where e/ = 1 + ^Q 1 dP/dEf\E f =o Substituting the num- 
bers, we find the theoretical value e/ — —525 for equilib- 
rium conditions while the experimental one found from 
Fig. lb (raw data for 2kHz, Ref.0) is e f = -680, i.e they 
are pretty close. According to Eq. JjJ, the absolute value 
of the negative slope for / = 9nm should be about two 
times larger than for I = 5nm, and one sees from the 
Figs. lb, 2 that this is indeed the case. If this agreement 
is not accidental, it indicates that in the films studied by 
Kim et al. the properties of the domain structure may 
be mainly defined by the electrostatics and not by some 
specific structural features of the samples, etc. The neg- 
ative slope of the hysteresis loops is apparently observed 
for the first time. 

To reveal specific features of domain structure in 
thinnest films we take into account that at small enough 
thickness the system should be in the paraelectric phase 
(as shown explicitly below) and consider loss of stability 
of this phase when the film thickness increases, i.e. a 
paraelectric-FE phase transition with thickness. The loss 
of stability is signalled by appearance of a non-trivial so- 
lution of the equations of state which can be either homo- 
geneous (SD) or inhomogeneous (MD). A homogeneous 
solution of 10) is possible at I > lh, where lh — L /e \A, 
at room temperature ^ T = 3.5nm. The domain instabil- 
ity means the appearance of a solution in form of a "po- 
larization wave' ' P = P k {z)e ikx [HEi of the (linearized) 
equation of state with the gradient term included: 

AP - gV\P = E, (8) 

where = d 2 + d 2 is "in-plane" Laplacian. In the case 
of metallic screening, this gives the following conditional : 



X t&n - xkd = e±k/e e \/k 2 + \ 2 , (9) 

where \ 2 — — e e±(A + gk 2 ) > 0, e± is the dielectric 
constant in the direction perpendicular to ferroelectric 



axis in the plane of the film. The case of interest to us 
is kX <C 1, easily met for metallic electrodes. We assume 
(and check validity later) that e±Xk/e e x ^ 1- Then, the 
equation simplifies to x/cd = 7r, the same as in FE film 
without electrodes or with a dead layer. We then find the 
maximal value (the highest temperature) of Ad = —2gk 2 
at k — k c where this equality is first met and domains 
begin to form: 

- A d = 2gk 2 = H/e l, k c = (ir 2 /e ± e a gl 2 ) 1/4 , (10) 

where £ = 2tt y/geo/e± is the characteristic length scale. 
Now, checking the assumption that we used to solve the 
Eq. J^J, we see that it boils down to Xe 1 / 2 '/e e (e .g) 1//2 > 1. 
Using values of A, e e from Ref. |lj, the value of g from 
Ref. [13, and calculating e^using the coefficients of Ref. 
Pj, we find that this condition is satisfied. From Eq. llOfl . 
we obtain the following critical thickness for domains at 
room temperature: 

If* = H/\A d \e Q ~ 3.1 nm. (11) 

Since ^ T < , the phase transition is into a MD state. 
The spatial distribution of a spontaneous polarization is 
near sinusoidal at I > ^ T . Higher harmonics develop 
with increasing thickness, and the polarization distri- 
bution tends to a conventional structure with narrow 
domain walls. But at small thicknesses it is nearly si- 
nusoidal, and one can expect weaker pinning compared 
to thicker films. It is hardly surprising that the em- 
pirical Merz's formula obtained for conventional domain 
structure does not apply to a sinusoidal one. The half- 
period of the sinusoidal domain structure can be esti- 
mated as o^ T ~ 7r/fc c = 1.7 nm at the transition and as 
a RT = 2.2nm for I = 5nm. 

It is instructive to consider the phase transition with 
thickness at zero Kelvin, where we get e± = 408, £ = 
0.08A, 4° K) = 1.6nm, lf K] = 2.5nm. The last result 
(homogeneous critical thickness of 2.5nm) is remarkable, 
since it practically coincides with the ab-initio calcula- 
tion for the critical thickness of 2.4nm in Ref. |14| . The 
ground state of the film is, however, not homogeneous but 
multidomain, and the domain ferroelectricity appears in 
films thicker than ^° K ^ = 1.6nm, which is the true "crit- 
ical size" for ferroelectricity in FE films in the present 
study. 
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